Statistics of extremal intensities for Gaussian interfaces 
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The extremal Fourier intensities are studied for stationary Edwards- Wilkinson- type, Gaussian, 
interfaces with power-law dispersion. We calculate the probability distribution of the maximal 
intensity and find that, generically, it does not coincide with the distribution of the integrated 
power spectrum (i. e. roughness of the surface), nor does it obey any of the known extreme statistics 
limit distributions. The Fisher- Tippett-Gumbel limit distribution is, however, recovered in three 
cases: (i) in the non-dispersive (white noise) limit, (ii) for high dimensions, and (iii) when only 
short- wavelength modes are kept. In the last two cases the limit distribution emerges in novel 
scenarios. 

PACS numbers: 05.40.-a, 05.70.-a, 68.35.Ct, 89.75.Da 



I. INTRODUCTION 

Extreme value statistics (EVS) has traditionally found 
applications in the analysis of environmental and engi- 
neering data, such as water level fluctuations or appliance 
lifetime Recently there has been a surge of activity 
from physicists in modeling phenomena that are governed 
by the extremal values of some random variable. Exam- 
ples range from the low-temperature state in spin glasses 
[3] , a variety of disordered systems with traveling fronts 
(see yl for a brief review), and modeling corrosive frac- 
ture Q , to depinning of surfaces (jj , to relaxation in 
granular materials |jj . 

A central notion in EVS is the existence of limit dis- 
tributions, that is, the extremal value in a batch of a 
large number of independent, identically distributed (i. i. 

d. ) random scalar variables obeys a probability density 
function (PDF), which is limited to one of three main 
types P, |3, 13 It is the tail of the PDF of the origi- 
nal variable, the parent PDF, that determines which one 
of these three classes the EVS will belong to. One dis- 
tinguishes the Fisher-Tippett-Gumbel (FTG) distribu- 
tion, for a decay faster than any power law, the Fisher- 
Tippett-Frechet (FTF), for power tail, and Weibull, for 
power law at a finite edge. The limit PDFs can even be 
cast into a single function with a parameter whose differ- 
ent ranges correspond to different traditional classes, see 

e. g. Ref. 0. 

The appearance of EVS limit distributions in physi- 
cal systems is in itself interesting and yields a new tool 
for the description of those systems, see e.g. On the 
other hand, there are instances when the quantities in 
question are either dependent, differently distributed, or, 
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have several components and thus generically do not be- 
longto any of the aforementioned EVS classes, see e.g. 
H UR 113 Closely related to our present subject are 
the studies on the distribution of local |l2|, [lj] and 
global extrema of the height of surfaces. Whereas 
there is no general mathematical theory for the EVS of 
non-i. i. d. quantities, quite a few special cases have been 
clarified, Refs. 0,0 has some examples of them. 

A possible connection of EVS to non-Gaussian fluc- 
tuations of spatially averaged, or, global, quantities in 
critical (stro ngly correlated) systems has recently been 
raised [MilpJlHjI^Jj^J^JglTl in relation to both exper- 
iments glllllllliiliiand simulations jHllllSl 
for many-body systems as well as for environmental 
data |2jj, [3(j- There are two striking characteristics of 
this ensemble of observations. Firstly, the form of the 
PDF is similar for these observables 0, |2(| . and sec- 
ondly, these functions bare a strong resemblance to the 
FTG distribution. There is some variation in form and 
there are some clear exceptions, but this approximate 
universality has been the subject of considerable debate 
in the literature. 

More concrete connections to EVS have been found in 
studies of the fluctuations of the roughness, or, width, of 
interfaces. The non-Gaussian nature of these fluctuations 
was first discussed in |3l| . and in 

families of PDFs were identified, all characterized by a 
single maximum, positive skewness, an exponential, or 
near exponential, tail for large fluctuations above the 
mean, and even faster decay below it. In Gaussian in- 
terface models relation to the FTG distribution has been 
observed in several instances. Namely, the magnetization 
distribution of the 2d XY model, in the low-temperature 
phase where vortices are rare, has been related to the 
roughness distribution in the stationary 2d Edwards- 
Wilkinson interface, and strong resemblance has been 
found to a formal generalization of the FTG function [26| . 
In another work on the 2d XY model, the FTG func- 
tion was numerically observed for extremal mode selec- 
tion [3^|. A strong connection to EVS comes from an 
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analytic result on a Id interface model with long-range 
interactions, corresponding to a Gaussian noise with 1// 
power spectrum, where the roughness of the interface 
turns out to be exactly of the FTG form [2l|. 

Given the manifold occurrence of the FTG shape, 
we are led to study the relation between EVS and the 
roughness of Gaussian interface models. These models 
can be considered as generalizations of the stationary d- 
dimensional Edwards- Wilkinson interfaces, specified by 
the exponent a in the power-law dispersion of Fourier 
modes (for a = 2 the Edwards- Wilkinson case is recov- 
ered). Due to their integrability, they are ideal starting 
point for the study of more complex, strongly fluctuat- 
ing systems. As the models diagonalize into statistically 
independent modes, and the roughness equals the sum of 
the Fourier intensities, it is straightforward to address the 
problem of possible connection to EVS. Namely, one can 
pose the question that originally motivated our research: 
do the modes with the largest intensities dominate the 
PDF of the roughness? Here we show that the answer to 
this question is no. In fact, this negative result, the first 
main conclusion of this paper, is immediately apparent 
once we construct the distribution of the maximal in- 
tensity an expression that turns out to be a generalized 
form of the Dedekind function. Besides being different 
from the distribution of roughness, this function generi- 
cally also deviates from the FTG or any known EVS limit 
distributions. The reason for the deviation from conven- 
tional EVS is that the independent modes have gapless 
dispersion, and are strongly non-identically distributed 
even for higher wavenumbers. As a result, for each real- 
ization of the interface, the largest intensity comes from 
one of only a few soft modes. The EVS is therefore ef- 
fectively coming from a finite-size system, even in the 
thermodynamic limit. 

The main body of the paper concerns the study of the 
maximal intensity PDF, which is found to depend on d 
and a. While there are no finite critical d c or a c , mark- 
ing thresholds to the known EVS limit distributions, we 
discover three limits where FTG statistics sets in: (i) 
a — ► 0, (ii) d — > oo, (iii) when only hard modes with wave 
vectors beyond a diverging radius R are kept. Common 
in these cases is that the thermodynamic limit can be 
taken beforehand, N — > oo, where N is the total number 
of modes, so one chooses the maximal intensity out of 
an infinite set, in contrast to the conventional EVS pro- 
cedure. The extremal value PDF becomes degenerate in 
all three cases in the sense that the standard deviation 
shrinks to zero on the scale of the mean. This is similar 
to the traditional FTG scenario when such degeneracy 
appears for TV — > oo. Now, however, in each case spe- 
cial scaling should be applied to resolve degeneracy and 
reveal the FTG distribution. It is suggested that the 
FTG limit in case (iii) is responsible for the numerically 
observed fit to the FTG function in Ref. |35|. 

Finally, we turn to the question how the EVS changes 
by considering a different choice of expansion functions 
(modes) . This study further illuminates the fact that the 



EVS is not generically described by any of the traditional 
EVS limit functions. Since the dominant contribution to 
the EVS comes from a few modes, the EVS will depend 
on the specific expansion functions. We find analytically 
a different family of extreme value PDFs for the inten- 
sities if the interface is expanded in sines and cosines, 
albeit the overall shape goes close to that of the PDF of 
the maximal Fourier intensities. 

We section the paper as follows. In Section [H] we de- 
fine the Gaussian model (|II A|) . present the basic formu- 
las for EVS for maximal Fourier intensities Ijll B|) , then in 
111 CI evaluate the EVS for the model showing 1// noise. 
General 1/ f a noise is considered in III Dl The case of 
arbitrary substrate dimensionality is treated in Sec. II I II 
with the general formulas for the extremal distributions 
calculated in ITTTAl In Sec. ITTTBl the EVS for the mag- 
netization modes for the XY system in the frame of the 
Gaussian model is investigated. The emergence of the 
FTG distribution in three limiting cases is described in 
IIVI The EVS of square amplitudes from the expansion 
in sines and cosines is discussed in [3 Sec. IVIl is devoted 
to a comparison to the distribution of the roughness, and 
Sec. IVIll contains the concluding remarks. We give the 
list of abbreviations in Appendix ^ the small- z asymp- 
totes of the extremal value PDFs are derived in Apps.lBl 
and IO and the finite- TV correction to the extremal dis- 
tribution is considered in App. [Dj Some details for the 
calculation of the PDF in the white noise limit is clari- 
fied in App. El ^d in App. [0 we determine, for a > d, 
the initial asymptote of the PDF of the roughness for 
comparison. 



II. GAUSSIAN SURFACES IN ONE 
DIMENSION 

A. 1// Q noise 

The interface at position a; on a one-dimensional (d = 
1) substrate is characterized by the height h(x). Replac- 
ing x by time t, h(t) can be thought of as the distance 
from the origin of a random walker. 

The roughness, or mean-square width, is given by 



w 2 {h) = [h{t)-h]\ (1) 

where over-bar denotes the average of walks, or interfaces 
on the interval < t < T 

F=±J%(t)dt. (2) 

The Fourier decomposition of the interface gives the am- 
plitudes of the fluctuating modes we are interested in 

JV 

h(t)= c n e 2mt ' T , c_ n = c*. (3) 

n=-N 
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Here h(t) is defined on Nq equidistant points (t = kAt, 
T = N At) and we introduced the notation N = (N — 
l)/2, with No assumed to be odd, before finally taking 
the thermodynamic limit N — > oo. The roughness is then 
the integrated power spectrum 



1 N 
w 2 {h) = — l c « 

71=1 



(4) 



In the models we consider, the time signal h(t) has pe- 
riodic boundary condition, it exhibits 1// Q power spec- 
trum, and the modes are independent Gaussian variables. 
The path probability is given by 



P[h(t)] oc exp[-S[h(t)]], 



(5) 



where the action S is given in terms of Fourier intensities 



|c„| 2 as 



JY 



5 = ( 7 T 1 - Q ^n c 



(6) 



n=l 



with exponent a defining the l/f a noise spectrum. The 
do is a parameter setting the effective surface tension 
and the power of T originates from dimensional consid- 
erations |34| . 

The PDF of the roughness for the Gaussian model has 
been extensively studied in [2(3, El 0- For a = 1 it 
was found analytically to coincide with the FTG function 
pi), one of the limit functions of EVS. Nevertheless, the 
roughness is not a priori an extremal quantity and it 
remains an open question why it obeys the EVS. One 
may argue that the sum is dominated by the softer modes 
whose extremal values possibly determine the cumulative 
behavior jlSl l35| . Motivated by this problem we study 
the largest contribution to the sum Q. 

An alternative choice of Fourier coefficients would be 
the set (lZec n ) 2 , (lmc n ) 2 , as studied in [33. These are 
the coefficients from the expansions by sines and cosines, 
and while this does not seem to be a physically significant 
modification, the EVS for this set is, in general, quanti- 
tatively different from that of the intensities |c„| 2 . Nev- 
ertheless, the PDFs for the two sets are similar in shape 
and lead to the same physical conclusions, so in order to 
keep unity of discussion we focus on the set |c„| 2 along 
most of the main text, and summarize the results on the 
(lZec n ) 2 , (Tmc n ) 2 at the end. 



B. Extreme value distributions 



We calculate the probability that the maximal inten- 
sity (cnl^^ is T a ~ 1 z /o"o, so we shall actually determine 
the scaling function of the EVS. We denote the PDF for 
that extreme value by P a (z), and the cumulative, or, in- 
tegrated probability distribution function (IPDF) by 



M a (z) = / P a (y)dy. 



(7) 



Since M a (z) is the probability that none of the \c n \ s 
exceed T a ~~ 1 z/a , we can express it as 

N f \c n \ 2 <z 

M a (z) = ]J A„e-" a|c " |2 dftec„dZmc„, (8) 

n=l Jo 

where A n — n a /tt is the normalization constant for the 
PDF of the nth mode. After evaluating the integrals in 
® we find 



N 



M a (z) = n 1 - 



(9) 



In fact, for a — I the product © is known from the 
defining formula of Dedekind's eta function, the latter 
also containing an extra power prefactor [3jj. Despite 
the difference in the prefactor, we will refer to JjJJ as 
a generalized Dedekind function. Differentiating 10 we 
arrive at the expression for the PDF of the maximal in- 
tensity 



A? 



(10) 



71=1 



The IPD F and the PDF vanish for z < 0, the above 
formulas I|9I10|I are understood for non- negative z. 

Three limits are easily understood here. Firstly, if a = 
0, i. e. the modes are identically distributed, we find 

M {z) = {l-e-*) N . (II) 
Then the change of variables z = ay + 7 + IniV, with 

a = v/C(2) = tt/V6, (12a) 
7= lim (Vn -1 -lnivV (12b) 



where £ is Riemann's zeta function and 7 Euler's con- 
stant, yields the FTG distribution in the N — > 00 limit 



M (z(y)) -> M FTG (y) = e~ e ' ay ~\ (13a) 
P {z(y))z\y) -> P FTG (y) =ae- a ^- e ~ aH ~ T (13b) 

The constants in (|12fl were used to scale the FTG distri- 
bution so that the mean becomes zero and the variance 
one |37| . Secondly and thirdly, in both limits a — > 00 
and z — > 00 only the mode n — 1 matters and one finds 



M a (z 
P a (z 



00 ) 
00) 



M^z) = 1 
Poc(z) = er 



(14a) 
(14b) 



Already at this point we can draw one of the main con- 
clusions of the present paper, namely, that for general a 
the extreme value PDFs are none of the known EVS limit 
functions. This immediately follows from the fact that 



4 



the negative-z and large-z behavior of (|9I10() is incom- 
patible with those of the limit functions for the statistics 
of maxima Q, Q ■ 

The breakdown of validity of the traditional EVS limit 
distributions can be ascribed to the fact that, for a > 0, 
due to the n a dispersion, the distributions of the individ- 
ual modes are sufficiently different. Here "sufficiently" 
needs to be emphasized, because one can conceive sets 
of different parent PDFs for the modes j^l that lead to 
say the FTG function. Such is the case of a dispersion 
that has a gap at the origin or goes to a constant for high 
frequencies. Conversely, one should not be surprised by 
the appearance of special PDFs different from the known 
limit distributions, because actually for any given PDF 
one can choose sets of different parent distributions yield- 
ing that PDF in EVS 0. 

It is easy to convince oneself that for a > in formulas 
(|9I10I) the limit N — ► oo can be taken and the resulting 
distribution has finite moments. Indeed, because of dis- 
persion the exponential PDFs of individual modes decay 
increasingly fast and so no singularities appear for large 
N. In conventional EVS, limit distributions arise be- 
cause the moments of variable z scale in singular ways 
in the thermodynamic limit, thus removing most of the 
details of the distribution of the original random vari- 
ables. Now, however, there is no such singular scale and 
the strong dependence on the statistics of the individual 
modes through the parameter a remains. This gives us 
a physical intuition about why the EVS for intensities of 
interfaces with dispersion is not described by any of the 
traditional limit distributions. 

Furthermore, we can also immediately assert the dif- 
ference between the PDF of the roughness (@J, studied 
in Ref . j^ij , and the PDF (|10|l of the maximal intensity 
component in An obvious deviation is in the physi- 
cal scales, namely, while the maximal intensity is always 
of order T 1_Q , the mean roughness is like that only for 
a > 1, diverges logarithmically for a = 1, and stays fi- 
nite for finite sampling times At if a < 1. Even when the 
scaling is the same, for finite a > 1, the two PDFs clearly 
differ (the generating function of the roughness PDF is 
given in |34|. whose Laplace transform is not Eq. l|l(J[l h 
and became the same only in the limit a — > oo, when the 
ii = 1 mode dominates. More detailed comparison will 
follow in Sec. EH 




FIG. 1: Distribution functions for the a = 1 case evaluated 
numerically according to Eq. IjlOfl for various TVs, each func- 
tion scaled to unit average. The convergence is apparently 
fast for practical purposes. 



0.01 



0.001 




FIG. 2: The same as Fig.Qon a semilog scale so as to better 
display the regions where -Pi (2) is small. 



be approximated as 

l/z 

M 1 (z)^l[(nz) 



-l/z 



(15) 



C. 1/f noise (a = 1) 

While for large z the PDF P\{z) is pure exponential as 
given in (|14b(l . for general arguments one should resort 
to numerical evaluation. The function l|lUfl is shown for a 
sequence of Ns in Figs.nand[3 For N > 9 the curves are 
hard to distinguish by the naked eye, the approximation 
by finite N thus converges fast for practical purposes. 

For small- z , the expected nonanalytic behavior can be 
estimated in the following way: the product in (|10|l can 



where we have used Stirling's approximation for the fac- 
torial. While here we took the terms with n > l/z as 
one, this gives us a first hint of the expected functional 
form. A more precise calculation in Appendix [5] shows 
that the large-n region gives a contribution of equal or- 
der, and in the end the true asymptote differs from i|15|) 
only in a scale factor in the exponent 



~ \l~ e"" 2/6 ", (16) 
see Eq. (|B14p with (|B12|) for a = 1 . Derivation by z gives 
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the small-;? asymptote of the PDF 

-/27T7T 2 



6z 5/2 



(17) 



As demonstrated in Fig.|3 the above expression is correct 
at small z and provides a reasonable approximation over 
most of the ascending part of P\{z). 



0.1 



0.01 



0.001 



0.0001 




FIG. 3: The PDF Pi(z) from Eq. (UJ with N 
small- z asymptote as given by Eq. 1171 . 



100 and its 



We conclude the case of the 1// noise by noting that, 
the PDF for the roughness scales logarithmically in 
N and approaches the FTG function |2l|. In contrast, as 
shown above, direct extremal selection of the constituent 
intensities in Q leads to a nonsingular PDF, related to 
the Dedekind function, even in the thermodynamic limit. 
In short, we have the FTG distribution for the rough- 
ness, which is not an extremal quantity, while the EVS 
is not describe d by the FTG function. Thus the question 
raised in Ref. |2l|, namely, what kind of extremal value 
selection can possibly be responsible for the FTG distri- 
bution of the roughness, has not been answered and is 
left to further explorations. 



the IPDF is given in Appendix IBI in (|B14|) . whence dif- 
ferentiation yields the PDF 



P a {z) 



(2n) a / 2 c(a) 

z 3/2+1/q 



■ exp 



z l/ a 



(18) 



where c(a) is given in Eq. I|B12(I . One can see that the 
above form becomes singular only in the a — > limit. 
While for any finite a > the asymptote is obviously 
incompatible with the FTG function (|13b|l . surprisingly, 
expression (I18fl corresponds to the generalized FTF func- 
tion for k th maximum Q in the special case when k, 
the FTF power parameter fj,, and a are related through 
fi = 2(fc— 1) = 1/a. The function l(T%|l does not equal the 
EVS distribution for larger z, it is not even normalized 
to one, so this coincidence does not contradict the claim 
that the extreme value PDFs P a (z) are none of the known 
limit distributions of EVS. The FTF class is not expected 
to be of relevance here anyhow, because the parent PDFs 
of the constituent modes decay exponentially. 

When one evaluates P a {z), the convergence in N is 
important. The iV-dependence has been determined in 
App-El whence it is apparent that convergence is fast for 
a of unit order and larger, but slows down for smaller a. 
We give here the asymptote in the region of slow decay, 
a < 1, 



N 

P a (z) - P N , a (z) » M a (z) — e 

a z 



(19) 



where we kept the N index for P when only N modes 
were counted. 

The results of the evaluation of P a (z) are given in 
Figs. 01 and [SJ If displayed on the scale when the mean 
is set to one, see Fig. 0J for vanishing a the PDFs de- 
velop a singularity as eventually they approach the Dirac 
delta. Slow convergence in N for small a is also demon- 
strated, already for a = 0.1 one has to go up to exceed- 
ingly large TVs to get a satisfactory approximation for the 
PDF. When the PDFs are shown with variance scaled to 
one, as in Fig. [31 the functions remain nonsingular and 
tend towards the FTG distribution. 



D. General a 

Here we go beyond the 1// spectrum and investigate 
the noise for general a. The special cases a = and 
a = 2 correspond to white noise and the Wiener process 
(ordinary random walk), respectively. The PDF of the 
roughness for Gaussian noise for periodic and "window" 
(bulk) boundary condition has been studied in [3(], 13 > 
and no a other than 1 has been found for which the PDF 
coincided with any of the known limit PDFs of extreme 
statistics. 

From the previous discussion it is clear that it is only 
for a — a c = that EVS is given by the FTG function 
in finite dimensions. Concerning the small-z asymptote, 



III. GENERAL DIMENSION 

A. Extremal intensities 

It is easy to generalize the above calculations to sur- 
faces h(r) defined on a hypercubic lattice substrate of 
dimension d and edge length L. We retain the periodic 
boundary condition for substrate, so the natural expan- 
sion of h(r) involves again Fourier modes. The probabil- 
ity of a surface (J5J is then characterized by the effective 
action [U E3 

JV 

S = a L d - a J2' H a K\ 2 - (20) 
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z/(z) 

FIG. 4: Distribution functions for d — 1 and various as. 
We rescaled z by its mean (z) to get a variable of unit av- 
erage. For a = 0.1 we show the sequence of approximants 
with N = 10 5 , 10 6 , 10 7 , peak moving upwards with increasing 
N, to demonstrate slow convergence. 



1.2 




-2-101234 



(z-(z))/Az 

FIG. 5: The PDF with the same as as in Fig. but z is 
rescaled now by its mean {z) and standard deviation Az to 
produce a variable of zero average and unit deviation. For 
a = 0.1 only the curve with the largest TV is shown. The peaks 
decrease with a. The PDFs tend to the FTG function (dashed 
line). The semilog inset magnifies the small-P regions. 



Here the rescaled wave vector n = (ni, ri2, rid) has 
integer components rii such that |nj| < N, furthermore, 
the mark ' implies that if an n is counted then — n is not, 
and the zero vector is excluded. (The halving of the Bril- 
louin zone is the consequence of the relation c n = c*_ n .) 
The above action for a — 2 corresponds to the stationary 
state of the original Edwards- Wilkinson model, while for 
a = 4 it gives the curvature-driven Mullins-Herring in- 
terface. For a general a = 2k, k integer, it is a Gaussian 



massless model with finite-range, and for a ^ 2k with 
interactions decaying like a power law. 

Denoting the maximal value of ao L d ~ a \c n \ 2 by z, and 
its IPDF by M a ^(z), each mode gives a multiplicative 
factor 1 — e~l™l z resulting in 

N 

M a , d (z) = ^'(l-e-l"!"*), (21) 

n 

whence the PDF is 

N ' \n\ a 

P a , d (z) = M a , d (z)J2 e |J h '_ r (22) 

n 

The IPDF H21(l can be considered as a further gen- 
eralization of Dedekind's original product formula [3g. 
As in one dimension, it is straightforward to show that 
for a > the above functions remain finite and involve 
finite mean and variance in the limit N — > oo. Thus 
again no singular scaling is necessary and so we do not 
expect that any of the known extreme value limit PDFs 
emerge for general d and a. However, again as in d = 1, 
for a — all independent modes become identically dis- 
tributed and we recover the FTG function after proper 
scaling by N. 

The large a limit of l)22[l. for any fixed dimension d, is 
determined by the contribution of the modes |n| = 1 

PooA^=de- z (l-e- z ) d ~ 1 . (23) 

For finite a where all modes count, we have closed forms 
only for the asymptotes. In the large-z limit, for any 
a > and fixed d, one obtains 

P a A z -> o°) ~ der z . (24) 

Since here, too, only the modes |n| = 1 matter, the large- 
z formula is obtained from the PDF IL'/il) for a — > oo, 
by taking z — > oo, which explains the fact that (|24|l is 
independent of a. For z — > we determined, in Appendix 
IB 2l the asymptote of the logarithm of the IPDF M a: d( z )i 
see Eqs. HB151IH16|I . One can easily convince oneself that 
the leading term for the logarithm of P a ,d{ z ) ^ s same 
as for M a ,d{z), so we have for small z 

For general z we have to evaluate expression (I22f) numer- 
ically. This poses no difficulties unless a w 0. 

B. The 2d Edwards- Wilkinson model (a = 2) 

One of the physically most relevant cases is d — 2 and 
a = 2, the Edwards- Wilkinson surface as originally de- 
fined. Since it is a massless Gaussian system, it also 
corresponds to the XY model in that part of the low 
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temperature phase where the effect of vortices is negli- 
gible |18|. There the magnetization distribution is pro- 
portional, with a change of sign, to the roughness of the 
Edwards- Wilkinson surface. Accordingly, the distribu- 
tion of the maximum amplitude of the magnetization 
fluctuations, after proper scaling, should be given by the 
PDF for d = a = 2. Fig. © shows the P 2) 2(z), to- 
gether with the FTG function normalized to the same 
mean and variance. While the fact that the two func- 
tions are different is obvious, here we demonstrate that 
the functions deviate significantly and in no range could 
they be mistaken for each other. The large difference 
here is of importance, because the PDF Pa,2(z) charac- 
terizes also the XY model, so this is an example when the 
FTG statistics is far from an EVS in a critical many-body 
system. 




z/(z) 



FIG. 6: Distribution function for the maximum Fourier in- 
tensity statistics for d = a = 2 and the FTG function scaled 
to unit mean and the same variance. 



the interesting region of z, where — In M(z) is not too 
large, z must be large, thus we can linearize the loga- 
rithms. Then we notice that for small a the terms in the 
sum change slowly with |n|, therefore we assume that the 
sum can be replaced by an integral. So we get 



\nM a ^ a ,d(z) 



2 ^ 

|n|>0 



7T 



d/2 



T(d/2) 



n d - l dn. (26) 



Asymptotic analysis of this expression is done in Ap- 
pendix [E| and results in 



where 



lnM Q ^o,dO) 



£ = -, 

a e 



■a(1 



d/a 



A = 



T(d/2) 



(27) 

(28a) 
(28b) 



valid in the region where z is £ in leading order. Intro- 
ducing the variable y by the linear transformation (for 
the constants a, 7 see Eq. 11121 ) 



z = t+-{ay + 1 + \nA), 

e 

and keeping y of order unity when a — > we get 

d/a 



A 



-ay-i 



(29) 



(30) 



up to terms vanishing with a. Since 1)27(1 was the leading 
term in the logarithm of the IPDF we get 



(31) 



IV. FTG LIMITING CASES 
A. Case (i): White noise limit a — > 

As we have seen, for a = the modes become identi- 
cally distributed in Eq. (|21|1 . thus, because of the expo- 
nential decay of the parent PDF for large z, one recovers 
the FTG function by proper scaling in N . One may, 
therefore, expect that when N = 00 is set first, the ex- 
treme value distribution converges to FTG for a — > 0, if 
proper scaling in a is applied. Such a convergence is not 
part of the standard theory of the FTG limit jlj,l3, so we 
show here how it comes about. A numerical demonstra- 
tion in d — 1 was given before on Fig. \5\ now we derive it 
analytically for any d, and determine the natural scaling 
of the variable z by a. 

Let us start out from the IPDF given by 1(21(1 . Ex- 
ponentiating the product into a sum, we realize that in 



which is the IPDF itTJal) for the FTG-distributed variable 
y. Since this IPDF has zero mean and unit variance, the 
linear transformation 1(291) is cquivalently 



yAz, 



(32) 



where (z) is the mean and Az the standard deviation 
of z up to terms vanishing for a — > 0. It thus follows 
that the scaled maximal intensities z have an average 
diverging proportionally to a" 1 , and they scatter in an 
O(l) region about the average. 



B. Case (ii): High dimensions 

We immediately recognize the FTG limit when d — > 00 
in the formula ((23(1 of the PDF for large a, Poo.d(z). In- 
deed, the PDF for the variable y — (z — 7 — logd)/a 
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then goes over to the FTG limit function (|13b|) . Below 
we study the d — ► oo limit for any fixed positive a and 
conclude that again FTG arises. This is far from obvi- 
ous a priori, because while in large dimensions there are 
many modes in each shell of constant \n\, so they have 
the same, exponential, parent distribution and, if we had 
only them, they would give rise to FTG, but on different 
shells the intensities are non-i. i. d. variables. 

Again we exponentiate the product (|21[1 into a sum, 
and notice that if d is large, the entropic weight quickly 
increases with |n|. Thus the dominant contribution 
comes from large |n|s, so we replace the sum by an inte- 
gral and again get formula (|26[) . which should be taken 
now for any fixed a > but in the limit d — > oo. Its 
asymptotic analysis is similar to that described for a — > 
in Appendix IeI and gives 



lnMo^ootV) 



1 



'2a 



d/o 



where 



i^T' 2 jl-q/2 



a e 



(33) 



(34) 



provided z is of the order of £. Next we introduce y 
through 



a ( In 2a 



Keeping y at order unity while d — ► oo, we get 



d/c 



-ay— 7 



V2a \z, 
whence (|33J) yields the IPDF 



(35) 



(36) 



(37) 



which coincides with the FTG distribution (|13a|l . 

We can thus conclude that in high dimensions the ex- 
tremal intensities belong to the FTG class. Furthermore, 
by comparing Eq. (|32[1 with (|35|) we can determine the 
average and standard deviation of the scaled maximal 
intensity z, namely, (z) oc d 1_Q / 2 and Az oc d~ a / 2 . In- 
terestingly, (z) diverges only for a < 2, while it shrinks 
to zero when a > 2, so the latter case is an example of a 
nonconventional FTG limit, when the maximal intensi- 
ties are small although the intensities themselves are not 
bounded from above. While it is plausible that strong 
enough dispersion can have the effect of an upper cutoff 
on the amplitude, the novelty here is the sharp transi- 
tion at a = 2, the only value when the characteristic 
maximal intensities remain finite and positive. For all 
a > 0, however, the scale of the standard deviation be- 
comes much smaller than that of the average, a feature 
of the conventional FTG scenario 0, . 



C. Case (iii): Hard modes 

We study the situation when the maximal amplitude 
is selected only from among those with |n| > R 3> 1, 
that is, the very hard modes. We consider arbitrary but 
fixed a and d. For we take the thermodynamic limit 
N — > oo, we are left with a single divergent parameter 
R. Whereas all hard mode intensities are quite small, 
they are of different scales and so are essentially non- 
i.i.d. variables. Thus special considerations are necessary 
to determine their EVS. 

The IPDF for the maximal hard mode amplitude is 



i(z,R)= n'(i 

\n\>R 



-|"| c 



(38) 



Obviously R — 1 gives the formerly studied IPDF l|21|) 
in the thermodynamic limit. The interesting region in 
z is where M(z) changes fast. At this stage we assume 
that in that region R a z is large, but we should check the 
result for consistency in the end. We can make from i|38|) 
a sum in the exponent, and since large |n|s are involved, 
we rewrite the sum into an integral. In leading order we 
obtain 

B R^~ a 

lnM atd (z,R -> oo) w exp (-R a z) , (39) 



where 



B = 



aT{d/2)' 



(40) 



Careful consideration of the compounding logarithmic 
singularities leads to the observation that in terms of the 
variable y, introduced by 

z = R- a Uy + 7 + ln— +dlnR-lnlnRj , (41) 



the expression (|3"9"|) is just 



-oj/-7 



up to terms vanishing 



with increasing R. That is, we have recovered the FTG 
function l|13a|) 



M aA (z, R —> oo) 



(42) 



where the linear transformation (|41|l is understood. 
Comparing (|41|l with i1M_?l) we find that the mean (z) has 
the leading singularity R~ a lni?, so in the relevant region 
of z the R a z indeed diverges, as assumed in the above 
derivation. 



D. Scales of singularity 

Summarizing the aforementioned limits, we found that 
FTG emerges (i) as a — > 0, (ii) for d — > oo, and (iii) 
when only hard modes R < \n\ with R — > oo are con- 
sidered. While for a = the intensities become i. i. 
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TABLE I: Order of the mean and standard deviation of 
the scaled maximal intensity z in the FTG limit cases (i-iii). 
For comparison the scales in the traditional scenario are also 
shown, see text for the parameters N, r. 



Case 






<^> 


Az = ij (z*) - (z) 2 


Az/ (z) 


(i) 


a - 


-» 


a" 1 


1 


a 


(ii) 


d- 


» oo 


d l- a /2 


d~ a ' 2 


d- 1 


(iii) 


R- 


-> oo 


R~ a lnR 


R~ a 


(lni?)" 1 


i.i.d. 


N - 


-* oo 


(\nN) 1/r 




(In AT) -1 



d. and so FTG should be expected, they are a priori 
non-i. i. d. in the cases (ii,iii). Nevertheless, for d — ► oo 
we found that a shell of practically i. i. d. modes in 
the Brillouin zone becomes dominant, so this essentially 
explains why one of the known EVS limit distributions 
emerged. On the other hand, in the case of hard modes 
we do not see i. i. d. intensities grouping, thus presently 
we lack an intuitive explanation for FTG. Remarkably, 
however, a common feature of all the above cases is that 
the distribution narrows down to a scale smaller than 
that of the average, a phenomenon also present in the 
traditional FTG scenario. In Tabled we summarize the 
scales of the mean and standard deviation of the max- 
imal amplitude, and compare them to the scales in the 
conventional limit when FTG emerges. In the latter case 
we consider batches of i. i. d. variables, whose parent 
IPDF approaches 1 like exp (— a z r ) for large z. The di- 
vergent parameter is then the number N of variables in a 
batch. Then the statistics of the maximal values within 
the batches becomes of FTG type ||, and straightfor- 
ward calculation yields the scales given in the last row of 
the table. 



imum of lZe 2 c n and Im 2 c„ as 

N , rV^ rVi 
M a . d (z) = TT / / 4e" |nncn| dftec^dZmcn. 

Then one straightforwardly gets 

N 

M a , d (z) = II erf 2 (v4^) , (44) 

n 

whence 

p a ,d = — i^^y^ — ?~j \ ■ ( 45 ) 

Vttz ^ erffyH^J 

This will be the basis for numerical evaluation where we 
go up to an N where the curve in the figure visibly sta- 
bilizes. 

For a — ► oo the \n\ = 1 modes numbering 2d dominate, 

so 

P^z) = M^erf 2 ^ 1 ^), (46) 

The large- z formula for any a is also determined by the 
softest modes, so it is just the asymptote of IpEri)) 

1dp~ z 

P a . d {z -> oo) « -^U, (47) 

while for z — > the logarithm of the PDF has the same 
leading term as the IPDF, given in Eqs. HC1IIC2II . 

For the sake of demonstration we consider a = d = 1, 
and denote the corresponding PDF by Pi. For small z, 
Eq. (IG3p of Appendix IU1 with a = 1 gives the asymptote 



V. SINE AND COSINE EXPANSION 

Based on physical intuition one may suspect that, in 
the absence of the traditional EVS limit for i. i. d. vari- 
ables, the extremal amplitude PDFs will depend on the 
choice of expansion functions. Given the periodic bound- 
ary conditions for the surface, another natural choice of 
expansion functions are the sines and cosines, whose co- 
efficients squared are (72.ec™) 2 , (Xmc n ) 2 . Below we show 
that indeed a new family of PDFs arise for the maximal 
square amplitude in this case, thus further illustrating 
deviation from the known EVS limit functions. More mo- 
tivation to look at these coefficients comes from the fact 
that the first EVS study in this area was done with such 
an expansion in Ref. |35j. There, the FTG distribution 
was found numerically when soft modes were discarded. 

When the sine and cosine modes are considered sepa- 
rately, we obtain the IPDF of z L a ~ d /<tq being the max- 



where c\ = 5(1, 1) = 1.3320405 was computed from l|C2|) . 
The full function and the asymptotic formulas are illus- 
trated in Fig. [7| Note that the deviation of the small- z 
asymptote from the exact PDF never exceeds 0.25 even 
for larger zs. We also displayed the PDF of extremal in- 
tensity Pi(z) from Fig. ^ which is given by a different 
formula and has a different mean (z) , but after rescaling 
goes surprisingly close to Pi- Nevertheless, the two func- 
tions can still be distinguished as shown on the inset of 
Fig. [7| This demonstrates that in the present case the 
EVS depends weakly on the expansion functions. 

Interestingly, in the three limits yielding the FTG dis- 
tribution for the intensities |ci| 2 in Sec. IIVI the FTG 
function is found also for (lZec n ) , (2mc n ) , following 
derivations similar to the ones in Sec. IIV CI So the FTG 
limits seem to be robust with respect to the choice of 
expansion functions. The FTG function found for the 
case of hard modes (the equivalent of the result in Sec. 
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and (|F9|) , that again can be written in the same form, so 
in Id both PDFs are asymptotically 



1.2 
1 





0.8 - 


< 


0.6 


3- 






0.4 




0.2 - 




- 




FIG. 7: The extreme value PDF Pi(z), computed from lHHt 
with a — d = 1 (full line), the asymptotes for small z from 
147H (dashed), and for large z I48H (short dashed). For com- 
parison the maximal intensity PDF -Pi (2) from Fig. is also 
shown (dotted); the two PDFs are very close, but distinct, as 
can be observed also near the maximum (see inset). 



II V Cf) explains the finding of |35( , where extremal square- 
coefficient statistics was studied numerically for the XY 
model, and the FTG function i|13b(l proved to be a good 
fit when soft modes were discarded from the batches of 
intensities before the selection of the maximal ones. 



VI. COMPARISON WITH THE DISTRIBUTION 
OF THE ROUGHNESS 

Interestingly, there is a resemblance between the scal- 
ing function of the extreme value PDFs for a > and 
those for the roughness [34| for a > d. Namely both 
types of PDFs have a single maximum, positive skew- 
ness, a nonanalytic initial asymptote, and a dominantly 
exponential decay (the leading term in the logarithm of 
the PDF is linear for large z). However similar they are 
qualitatively, the functions are different. 

A straightforward way to make the comparison quan- 
titative is to calculate the asymptotes of the PDFs. For 
small x the roughness PDF behaves as given by (|F16(1 of 
Appendix|Fl and the EVS PDF goes like It is inter- 
esting to note that for the roughness the critical a c = d, 
where the variance shrinks to zero on the scale of the 
mean, while for the EVS the same type of criticality is 
observed for a c = 0. Then both asymptotes can be cast 
in the common form 



lnP(y) cx y- d /( Q - 



(49) 



where for the roughness and for EVS one should under- 
stand y as z and x, respectively. 

Furthermore, in Id one can calculate the power pref- 
actor in front of the nonanalytic exponential, c. f. i|18|) 



P(y) « C iy ~ 



x exp 



C 2 



y 



1/(0- 



(50) 



albeit the proportionality constants in the two asymp- 
totes are different. What is more, in Id both PDFs be- 
come the FTG function in the limit a — > a c . Note how- 
ever, that the two kinds of PDFs have different asymp- 
totes for larger arguments, when a > a c . 

We can put the threshold behavior in a short form. 
Firstly, for both kinds of distributions there is a (lower) 
critical value (d/a) c where, for increasing d/a, on the 
scale of the mean the distributions become the Dirac 
delta, and this threshold is 1 for the roughness and 00 for 
the EVS. Then there is an upper critical value (d/a)™, 
that is the threshold for the respective classical limit dis- 
tribution, Gaussian for the roughness and FTG for the 
maximal intensity. For the roughness we have seen |34[ 
that (d/a)" = 2 and for the maximal intensity it is again 
(d/a)" = 00. So there is a region, 2 < d/a < 00, where 
the roughness is Gaussian, but the EVS is still not given 
by any of the known limit distributions of EVS, rather 
by the generalized Dedekind function. 

There is a significant difference also in the finite-size 
correction to the PDF. The PDF for the extremal ampli- 
tude converges essentially exponentially fast (see App.IdI 
for the Id correction formula), while the correction to the 
roughness PDF can be shown to be algebraic in N. 

In sum, there is a strong qualitative resemblance be- 
tween the shape of the roughness PDF for a > d and the 
EVS, furthermore, their initial asymptotes have similar 
functional forms, they are, however, distinctly different 
functions. It remains to be clarified whether the similar- 
ity has some deeper reason, or it is simply a mathematical 
coincidence. 



VII. CONCLUSION 

The results presented above bring us to a very definite 
conclusion concerning the relation between interface fluc- 
tuations and extremal statistics: the roughness PDF is 
not given by the largest mode in the Gaussian interface 
model. This is true, even in the case a = d = 1, where the 
roughness PDF is the FTG distribution [2]]], one of the 
known limit functions of extreme statistics. We showed, 
further, that the PDF was none of the known EVS limit 
distributions and depended on the statistics for the indi- 
vidual elements of the model. In addition, one expects 
the boundary conditions to influence the shape of the 
PDFs. It should be added that, to the extent that the 
Gaussian model is a universal family of massless models 
from the viewpoint of critical phenomena, the general- 
ized Dedekind PDFs are equally universal, depending on 
the dimensionality d and the dispersion parameter a. 

It is worth pointing out that there is an analogy be- 
tween the Gaussian distribution arising from the central 
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limit theorem, which applies for sums of random vari- 
ables with finite moments, and the extreme value limit 
distributions such as FTG, which is about the maximal 
of those variables. Both limit distributions are related to 
a large ensemble of independent, identically distributed 
objects. The appearance of a non-Gaussian PDF for the 
integrated power spectrum ( i.e. rou ghn ess), in Gaussian 
systems [lj| [lTJ |2j, |23|, |2(| lU HH, Liil Lid , is a conse- 
quence of strong dispersion, whence follows the strongly 
non-identical distribution of the modes. What we have 
illustrated in this paper is that it is this same dispersion 
that generically excludes the known limit distributions of 
extreme statistics. To refine the picture, we found that 
if Gaussian central limit statistics is excluded for the in- 
tegrated variable because of dispersion, FTG extremal 
statistics is also explicitly excluded, but the reverse is 
not necessarily true: there is the region of finite d/a > 2, 
where the integrated power is Gaussian but the maxi- 
mal Fourier intensity does not follow FTG. However, we 
found three border cases where, by singular scaling, the 
FTG distribution arises. In the limit a — > the am- 
plitudes obviously become identically distributed. For 
the other cases, that is, in the limit d — > oo and when a 
large number of modes are omitted near the center of the 
Brillouin zone, the number of contributing modes diverge 
but the dispersion across the zone remains. This result 
appears to broaden the scope of validity of the FTG dis- 
tribution in EVS to non-i. i. d. variables. Furthermore, 
the special scalings that led to the emergence of FTG 
in the border cases demonstrate new ways of extracting 
the known limit distribution even in the absence of the 
conventional FTG scenario. 

In summary, the extreme value PDF for Gaussian in- 
terface models are, for fixed choice of variables, a two- 
parameter family of functions with a nonanalytic part 
for small values, a single maximum, and an essentially 
exponential tail. The functions look qualitatively simi- 
lar, but vary quantitatively over the range of parameters 
d and a studied. This family of curves is different from 
that found for the roughness distribution for the same pa- 
rameters a,d [34(. We conclude that it is not possible to 
make a direct link between non-Gaussian roughness fluc- 
tuations and extreme values in these models. Given the 
multitude of recent observations of non-Gaussian fluctu- 
ations in experimental and in model strongly correlated 
systems, as well as in the use of interface models as phe- 
nomenological tools in describing such fluctuations, this 
seems like an important result. However, the relevance of 
extreme values in more complex, non-Gaussian systems 
remains an open question and it would be interesting to 
follow up this work with studies in this direction. 
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APPENDIX A: ABBREVIATIONS 

Below we summarize the abbreviations used through- 
out the paper. 

extreme value statistics EVS 

independent, identically distributed i. i. d. 

probability density function PDF 

integrated probability distribution function IPDF 

Fisher-Tippett-Gumbel FTG 

Fisher-Tippett-Frechet FTF 



APPENDIX B: SMALL- 2: ASYMPTOTE OF THE 
EXTREME VALUE DISTRIBUTION FOR THE 
INTENSITIES 

1. d=l 

Our starting formula is Eq. @ where, for small z, we 
expand terms with not too large n to leading order in 
n a z. Terms with larger ns must be considered without 
expansion. However, exponentiating them the product 
becomes a sum of logarithms, which we can replace by an 
integral as the terms vary slowly with n. Carefully treat- 
ing various corrections allows us to obtain the asymptote 
for ©. 

Let us separate the product for the IPDF, while 

z — > 0, as 



M a {z) = CxD, 



(Bl) 



where 
C 



n a z = z n " (n Q !) Q ,(B2) 



D = 



n (i-e--')«n 

71=1 71=1 

exp ( jr ln(l-exp(-n"z)) ) , (B3) 



\n—n a -\-l 



with (n a ) a z small but n a large, i. e. 



» n a > 1, 



A/a 



(B4) 



The first inequality allows the linearization of the expo- 
nential in i|B2|> , whereas the second one will enable us to 
use the Stirling formula in ljB2fl and replace the sum in 
(IB3(I by an integral. 

A short detour is necessary to see under what condition 
the correction to the term linear in z in each factor of 
(|B2|) can be neglected. Including the next term in the 
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expansion we have 



C w j| n Q z f 1 - \n a z\ 

n=l ^ ' 



«p E» 8 n 

\ n =i / n =l 



n z. 



(B5) 



The exponent in the prefactor goes like zn^ +1 . So the 
prefactor can be taken as unity if we use an n a that 
satisfies 



1 



7TJ(a+l) » Ha » X ' 



(B6) 



a condition stricter than (|B4|) . Then the approximation 
in 1|B2(I indeed gives the asymptote of C. 

Now we can calculate C from the r. h. s. of (|B2I) by 
using Stirling's formula 



(77 \ an a 
^) (27rn Q ) ( 



/2 



(B7) 



To calculate D we first estimate the error incurred 
when the sum is written as an integral: given a func- 
tion f(z), the terms in the sum of f(n) from n a + 1 to 
infinity can be written to leading order 



71+1 



1 



/(n+l)«j <Lr/(aO + -/ / (n + l). (B8) 

Performing the summation, we can replace the sum of the 
f'(n + l)s by an integral in leading order and we finally 
find up to the next-to-leading order 

£ /(» + !)«/ dx/( a :)--/K). (B9) 

„=„_ Jn Q 



This is the exponent of (|B3fl with 

f(x) = In (1 - exp(-x a z)) (BIO) 
whose integral can be written using u — xz x l a as 
1 



dx f( x ) = I du\n(l~ exp(-u a )) 

1 



A/a 



du In (1 — exp(— /)) 
du lnu. 



(Bll) 



In the last line we used the smallness of u: one can easily 
convince oneself, from (|B6J| that the higher order terms 
vanish. The definite integral is 



c (a 



du In (1 — exp(— u a )) 



(B12) 



and so from l|B9|) we get 

InD » + «n Q - L + In ((n Q ) Q z) . (B13) 

Finally, (JbTJ) and (|B13)l give 

«.w-2^-p(-^). (Bi4, 

This is the sought after asymptote for small z in the sense 
that its relative error vanishes for z — > 0. 

Remarkably, when we neglect the large n contribution, 
we find the correct l/z 1 /" dependence in the argument 
of the exponential in M a (z). However, the coefficient 
c(a) is only found by keeping terms with large n. Note, 
furthermore, that the inequality ljB6|) needs a > 0. Thus 
we should not be surprised that the asymptote for a — > 
does not relate to the tail of the FTG function (|13b(l for 
large negative argument. 

In summary, the leading term of In M(z) was produced 
by exponentiating the product and rewriting the sum as 
an integral. We shall follow this recipe below for general 
dimensions. 



2. Arbitrary dimension 

We consider here the asymptotes of the IPDF 12111 . 
The leading term of the logarithm of the IPDF M a ^{z), 
in the small-z limit, can be determined by transforming 
the sum over the modes of the Brillouin zone into an 
integral, because for z — > the terms in the sum change 
slowly. Thus we obtain, to leading order in z 



\n.M a , d (z) 



c(a, d) 

z d/a ' 



(B15) 



where 



,d) = -~f d d uln(l-cxp(-|ur)) 



dT 



(B16) 



For d — 1 we indeed recover the c(a) of (|B12(I . Note that 
(IB 15(1 does not give the full asymptote for M, it is only 
the leading term in the exponent. 

In order to estimate the next-to-leading term for 
lnM Q ( ;(z) one can calculate the correction arising when 
the sum is transformed into an integral. This is of the 
order 



In ; 



z (d-l)/a 



(B17) 



This correction diverges algebraically for d > 1, giving 
rise to a further exponential singularity in the IPDF. The 
asymptote of M is therefore not as simple as in the Id 
case (IB 141). 
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APPENDIX C: SMALL-z ASYMPTOTE IN THE 
CASE OF THE SINE AND COSINE EXPANSION 

In this Appendix we just summarize the results for the 
PDF of the maximal square-coefficient for the sine and 
cosine expansion functions. The leading term for the 
logarithm of the IPDF l|44|l comes from our making the 
sum an integral 



In M a , d (z) 



c(a, d) 
z d / a : 



where 



c(a, d) 



d d u In 



(erf(M Q / 2 )) 



(CI) 



(C2) 



In d = 1, after a calculation along a line similar to the one 
followed in Appendix IB II we obtain the full asymptote 



- , , (27r) Q / 2 V^ / c(a,l) 
M a ,i(z) w 1 ' n ^ exp ' 



2^ 



yl/c 



(C3) 



Note that the functional form is similar to the asymptote 
of the IPDF for the maximal intensity (|B14|t . but the 
constants are different. 



APPENDIX D: FINITE- AT CORRECTION 



If a > 1 then 



I N ,a,l(z) = J2 e ""° Z= E e ~"" Q ^ 



n=N+l 
oo 

< E « 

n=Ar+l 



n=N+l 



-n{N+l)"- 1 z _ £ 



-(JV+l) a z 



1 - e -(^+i) c 



(D5) 



The last formula, an upper bound for the asymptote, is 
just the first term in the sum /. The sum has positive 
terms, whence it follows that the first term is at the same 
time the asymptote of / for a > 1. Finally, in the case 
a < 1 we can rewrite the sum for / into an integral, 
because if in Eq. I|B8|) we substitute e~" z for f(n) and 
iV+1 for n a , the correction to the integral is negligible for 
a > 1 in the large N limit. The integral is approximately 
N l ~ a e~ N joLZ, this is then the sought asymptote for /. 
In conclusion, 

ie - ( iV+i)« Z) if a >l, 
I N , aA (z) { = e ~ Nz (e 2 - l)" 1 , if a = 1, (D6) 



N 



l-a„-N° 



' J 'a z, if a < 1. 



We can summarize the above results for various as such 
that if the summand in / decays slowly, one can replace 
the sum by an integral, and if it decays fast then the sum 
asymptotically equals its first term. 



Here we calculate the correction of the extreme inten- 
sity distribution for large N. Denoting now the finite 
product H21fl by Mjv,a,d(^) an d assuming TV to be large 
we obtain 



M Ni a. d (z) 



M a4 



|n|>JV 



|n|>JV 



,-i«r 



1-I^„i(4 (Dl) 



Hence 



Mjv )et) d(z) w M a ,d{z){X + lN,a,d{z)) 

PN,a,d( z ) ~ Pa,d( z ) + M a ,d{z) l'N, a ,, 



iz), 



(D2) 
(D3) 



where the second line was obtained by differentiation and 
we used the property that / is negligible next to I' for 
large N. From I|D1(I we surmise that the convergence 
to the N — oo functions is in essence exponentially fast 
in N a ; to be specific, we give below for Id the precise 
asymptote. 

Considering in d = 1 the a = 1 case, we have a geo- 
metric series as 



n=N+l 



-Ns 



(e z -iy\ (D4) 



APPENDIX E: ON THE FTG LIMITS 

Below we derive Eq. (|28|l from (|26|l . We shall replace 
the lower limit of integration 1 by n$ to show irrelevance 
of the precise setting of the lower limit. Denoting the 
integral in (|26J) by I and changing the integration variable 
to v = n a z we get 



1 



a z 



d/a 



e -v v d/ a -l dv _ 



(El) 



Since a — > 0, the lower integration limit becomes z, in- 
deed independent of no- The saddle point method gives 



7T 

2ad 



d 

zae 



d/a 



erfc 



(E2) 



When z « d/ae, the argument of the erfc becomes a 
large negative number, where the erfc is approximately 
2. Thus, including the prefactor before / in (|26|l . we 
recover (|28H . 

We have here an opportunity to test whether the sum 
in (|26|l was justly rewritten into an integral. We do 
not have rigorous results, but we know that the inte- 
gral must not be smaller than a single summand term, if 
it is, the integral representation cannot be accepted. A 



characteristic summand term now is e 



where 
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v = (d — I) /a is the saddle point in (|E1|) . so the sum- 
mand term is small. Since we consider zs for which the 
full integral for — InM is of 0(1), the integral is indeed 
much larger then a single characteristic summand term 
e~ v . This is a generic test that the integral represen- 
tation of a sum should pass, and we performed it in all 
pertinent cases in the paper. 



APPENDIX F: SMALL-z ASYMPTOTE OF THE 
ROUGHNESS DISTRIBUTION 

1. One dimension, a > 1 

We derive here the small- a; asymptote of the PDF of 
the roughness for the Gaussian model of 1// Q noise with 
periodic boundary condition. We consider the case a > 
1, when all cumulants of the roughness are of the same 
order. In the case of the Wiener process, a = 2, the PDF 
has a nonanalytic asymptote, which has been calculated 
in Ref. j^l). Here we find that nonanalyticity prevails 
for all a > 1. We begin with the simplest form, free of 
normalizing constants, of the generating function [34j 



G a {a) = n f 1 



71=1 



S 

n a 



(Fl) 



The average is given by —G' a (0) — ((a), where £ is Rie- 
mann's zeta function, so in order to obtain a PDF nor- 
malized to unit average, in the end we should rescale by 

C(«). 

The large-s asymptote of the generating function will 
give us the initial asymptote of the PDF. The calculation 
goes along the lines of Appendix [BJ For large s we can 
factorize G a as 



G a (s) = Ex F, 



(F2) 



where 
E 



n=i 



F = exp 



ih^) -riT =s " n " Ma ' (F3) 

71— 1 

OO \ 

E ( F4 ) 



n— n a +1 



where 



a V(«+i) 



> n a > 1. 



(F5) 



The first inequality ensures that (|F3|1 indeed gives the 
asymptote for E, see App. IBIfor an analogous estimate, 
and the largeness of n a enables us to approximate the 
sum in l|F4(l by an integral. We also calculate in F the 
nonvanishing corrections to the integral, again in a way 
similar to what was done in App. [BJ Then using the 
Stirling formula in l|F3(l for E, substituting the E and F 



into (F2j, and collecting all stray terms we wind up for 
large s with 

G a (s) « (2n) a/2 ^ exp (-s 1 /" g(a)) , (F6) 



where 



9{ot) 



/"CO 

/ du In (1 
Jo 



sin I 



(F7) 



The PDF of the roughness, normalized to unit mean, is 
obtained by inverse Laplace transformation 



Pa(x) 



ds 



G 



(F8) 



For small x the large-real-s region dominates, where the 
saddle point method can be used. It suffices to compute 
the saddle point from the exponential of jF6j), and then 
substitute its value into the *Js prefactor. Taking into ac- 
count also the quadratic deviation from the saddle point 
in the exponent, we finally obtain the small- a; asymptote 



exp 



R{a) 



x 



l/(a-l) 



where 



Q(a) 
R(a) 



(27r) a 3-fl(a) 



y/a-1 (aCCa)) 2 ^- 1 ' 
(a-l)s(a)^ 



a > 



C(a)< 



(F9) 

(F10) 
(Fll) 



2. Arbitrary dimension, a > d 

In general dimensions we can give here only the leading 
exponential for the initial asymptote of the PDF of the 
roughness for periodic boundary condition. We consider 
the case a > d, else there is no sense in speaking about 
the asymptote for small roughness on the scale L d ~ a , see 
Ref. [341. T he generating function yielding unit average 
is from |34| 



G a ,d(s) 



n 

\n\>0 



C(a,d) \n\ 



-1/2 



where 



t(a,d) = l £ 



n 



(F12) 



(F13) 



|n|>0 



is a (i-dimensional generalization of the Zeta function. 
Exponentiation of the product to a sum and transforma- 
tion of the sum to an integral gives the leading exponen- 
tial term 



In G at d(s) fa - 



C(a,d) 



d/a 



g{a,d), 



(F14) 
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where 



with 



g{a,d) 



d d u ln(T 



T d/2+l 



dT(§) sin(^) 



(F15) 



R(a,d) 



Performing the inverse Laplace transformation by the 
saddle point method, one arrives at the small- a; asymp- 
tote of the PDF as 



T d/2+l 



«r(|) sin(^) a^dy/c 



,(F17) 



lnP ajd (a;) 



R(a, d) 



(F16) 
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